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ABSTRACT 



We present a new near-infrared imaging survey in the four CFHTLS deep fields: the WIRCam Deep Survey or "WIRDS". WIRDS 
comprises extremely deep, high quality (FWHM ~ 0.6") J, H and A* s imaging covering a total effective area of 2.1 deg 2 and reaching 
AB 50% completeness limits of « 24.5. We combine our images with the CFHTLS to create a unique eight-band ugrizJHK s photo- 
metric catalogues in the four CFHTLS deep fields; these four separate fields allow us to make a robust estimate of the effect of cosmic 
variance for all our measurements. We use these catalogues in combination with w 9, 800 spectroscopic redshifts to estimate precise 
photometric redshifts (<x A - /(1+ -) < 0.03 at i < 25), galaxy types, star-formation rates and stellar masses for a unique sample of as 1.8 
million galaxies. Our JHK S number counts are consistent with previous studies. We apply the "BzK" selection to our gzK filter set 
and find that the star forming BzK selection successfully selects 76% of star-forming galaxies in the redshift range 1.4 < z < 2.5 in 
our photometric catalogue, based on our photometric redshift measurement. Similarly the passive BzK selection returns 52% of the 
passive 1.4 < z < 2.5 population identified in the photometric catalogue. We present the mass functions of the total galaxy population 
as a function of redshift up to z = 2 and present fits using double Schechter functions. A mass-dependent evolution of the mass 
function is seen with the numbers of galaxies with masses of M < 10 10 75 still evolving at z < 1, but galaxies of higher mass reaching 
their present day numbers by z ~ 0.8 - 1. This is consistent with the present picture of downsizing in galaxy evolution. We compare 
our results with the predictions of the GALFORM semi-analytical galaxy formation model and find that the simulations provide a 
relatively successful fit to the observed mass functions at intermediate masses (i.e. 10 < log(M/M Q ) < 11). However, as is common 
with semi-analytical predictions of the mass function, the GALFORM results under-predict the mass function at low masses (i.e. 
log(M/M G ) < 10), whilst the fit as a whole degrades beyond redshifts of z ~ 1.2. All photometric catalogues and images are made 
publicly available from TERAPIX and CADC. 

Key words, cosmology: observations — cosmology: dark matter — galaxies: large scale structure of the universe — galaxies: surveys 
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1. Introduction 

Today, observational evidence concurs that star formation in the 
Universe reached its peak within the redshift range 1 < z < 2, 
whilst ~ 50% - 70 % of mass assembly took place in the re dshift 



2003; 



range 1 < z < 3 (Connolly et aljfl 9971; [Dickinson etal. 
lArnouts et all l2007t iPozzetti et all 120071; iNoeske et al. 
In addition, observations have revealed the existence of a sub- 
stantial population of very massive galaxies at z > 1 (e.g. 



2007). 
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* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de l'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is based in part on data products 
produced at TERAPIX and the Canadian Astronomy Data Centre as 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- 
rative project of NRC and CNRS. 



Mulli s et alJl2005b IStanford et al.ll2005l l2006t IStott et aLlEoioL 
iFassbender et al.ll201 ll) . Such observations have acted to stimu- 
late significant progress in our understanding of galaxy forma- 
tion processes, in particular highlighting the importance of feed- 
back processes in shaping the pace and form of the growth of 
galaxies. Semi-analytical simulations have shown that the ad- 
dition of AGN feedback allow such models to predict both the 
substantial population of massive galaxies at z > 1 and the rise in 
the star-formation density to z ~ 2, whilst simultaneously match- 
ing the luminosity functi on and colour distribution of galaxies in 
the local Universe (e . g. iBower et al.l 120061; ICroton et al.ll 20Qi 
| De Lucia etal] l2006t jPe Lucia & Blaizotl 120071; iGuo & Whit* 
2008: iGuo et al.ll2011h . 



Although models (and hence our understanding of galaxy 
evolution) are continually improving and have had many suc- 
cesses, there are still limitations and inconsistencies. For ex- 
ample, semi-analytic models have difficulty simultaneously pro- 
ducing the total mass function of the massive galaxy popula- 
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tion at z > 1 and stellar growth in highl y star-forming sub- 
millimetre galaxie s observed at z ~ 1 (e.g. Baugh et al. 2005; 
iBower et alj|2006l) . In addition, there remain i nstances of obser- 
vation s challenging the model predictions (e.g . lMcCracken et al] 
120 lOt iGuo etalj|20Tll iHenriques etal]|2012l) . which highlight 
that our understanding of galaxy formation remains incomplete. 
Given the central role that intermediate redshift populations play 
in our understanding of galaxy formation and evolution as a 
whole, gathering deeper wide-area observations at this redshift 
range is of crucial importance. 

Below z ~ 1, the galaxy population can be relatively eas- 
ily identified via the optical features such as the 4000A break 
and a variety of absorption and emission lines. Above z ~ 
2, the ultraviolet spectral features are redshifted into the op- 
tical and the galaxy populat ion can be readily i dentified via 
the Lyman Break (e.g. [Steidel et all 1 19961 120031: ICooke et alJ 
2006; Bielbv et al. 201 lj ) and Lya emission (e.g. ICowie & Hu 
199a lOuchi et al.M2003b iGawiser et al.ll2006l) . However, to fol- 



low the galaxy population into the redshift range 1 < z < 2.5, 
near-infrared observations are essential as all spectral features 
move out of visible bands. Furthermore, the role of environment 
and gal axy formation process a t these redshifts is largely unex- 
plored dRenzini & D addi 2009). Near-infrared galaxy samples 
offer several well-known adv antages compared to purely optical 
selections (see, for example ICowie et ail (U994)); at z ~ 1 near- 
infrared selected galaxies are seen in the rest-frame optical. This 
corresponds more closely to a stellar-mass-selected sample and 
are therefore less prone to the uncertain effects of dust extinc- 
tion present in these redshifts for rest-frame optically selected 
samples. As the K s band k— corrections are insensitive to galaxy 
type over a wide redshift range, these samples provide a fairly 
unbiased census of galaxy populations at high redshifts (provid- 
ing that the extinction is not too high, as in the case of some 
submillimeter galaxies). Such samples represent the ideal input 
catalogues from which to extract targets for spectroscopic sur- 
veys as well as for determining accurate photometric redshifts 
an d making compariso ns with models. 

'■ et alJ (1 19961) carried out one of the first extremely 
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deep, complete K selected spectroscopic surveys and showed 
that star-forming galaxies at low redshifts have smaller masses 
than actively star-forming galaxies at z ~ 1, a phenomenon 
known as "downsizing". Stated another way, the sites of star- 
formation "migrate" from higher-mass systems at high redshift 
to lower-mass systems at lower redshifts; this anti-hierarchical 
behaviour seemed at odds with the hierarchical picture of galaxy 
formation and helped to stimulate refinement of theoretical mod- 
els of galaxy formation t o account for these observations. More 
recently, the K20 survey (Cimatti et al. 2002) reaching K s «s 21.8 
and the GDDS survey (I Abraham et al.l l2004) reaching K » 22.4 
provided further observations to test the existing knowledge and 
understanding of the galaxy formation process. The areas cov- 
ered by these surveys was small, comprising only ~ 55 arcmin 2 
and ~ 30 arcmin 2 in K20 and GDDS respectively. Improvements 
in recent years have been provided by the Ultra Deep Survey 
(UDS) segment of the United Kingdom Infr ared Telescope 
(UKIRT) Infrared Deep Sky Survey (UKIDSS: [Lawrence et all 
120071: iLaneetaLl 120071) . which covers an area of 0.77 deg 2 to 
a depth of K ~ 25 (AB) and the Cosmic Evolution Survey 
(COSMOS; IScoville et al.ll2007t iKoekemoer et~afll2007l) which 
covers a total area of 2 deg 2 to a depth of K — 24 (AB). Each 
of these surveys covers just a single field however, and cosmic 
variance effects are particularly important, especially for the red 
galaxy populations that are dominated by highly-clustered ob- 
jects. 



A number of studies have analysed the evolution of the build- 
up of stellar mass in galaxies via the stellar mass function, find- 
ing that the stellar mass density approximately doubles between 
z ~ 1 and z ~ (e.g iBundv et al 112 005. 20061 Borch e"tai]l2006t 
iPannella et alJ^OOCT . IIlbert et alJ (l2010h performed similar anal- 
yses out to redshifts of z ~ 2 and found that, for quiescent galax- 
ies, z ~ 1 represents an epoch of transition in their stellar mass 
assembly, showing that the stellar mass density for the quiescent 
population increases by over an order of magnitude from z ~ 2 to 
z ~ 0.8, but only increases by a factor of 2 between z ~ 0.8 and 
the present day. Combining this result with the morphologies of 



the 



< 



0.8 quiescent sample, they concluded that a dominant 



mechanism links the shutdown of star formation and the acqui- 
sition of an elliptical morphology in massive galaxies. Although 
these studies have provided significant insights into the evolu- 
tion of the galaxy population, building on these results requires 
further wide field deep NIR imaging to probe a fuller range of 
the stellar mass content across cosmic time. 

In this paper, the first of a series, we present such a sur- 
vey: the WIRCam Deep Survey or "WIRDS". We describe the 
field layout, observing strategy, data processing techniques and 
present a comprehensive quality assessment of the final released 



studies: Bielbv et al. 


2010; Cassataet al. 2011; Andreon 
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rCucciati et al.ll2012; 


Raichoor & Andreon 2012tlWhitaker et alJ 



2011; Willottet al. 2012. Both images and catalogues are pub- 
licly available from CADCQ and TERAPIX0 data centers. In 
our second paper, we will describe clustering measurements for 
mass-selected samples derived using these catalogues. 

This paper is organised as follows: Section [2] presents the 
survey strategy, observations, data reductions and quality as- 
sessments; Section|3]we present galaxy number counts and sub- 
classes of objects selected using various selection criteria in joint 
CFHTLS-WIRDS dataset, including BzK selections and separa- 
tion by type based on rest-frame UV colours. In Section [4] we 
present an analysis of the mass completeness of the data and 
present precise mass functions measured to a redshift of z — 2. 
Section[5]provides a summary and our conclusions. The WIRDS 
survey is designed to reach AB magnitudes of 24 in all near- 
infrared bands. The science rationale for this was the following: 
firstly, for the purposes of deriving accurate photometric red- 
shifts below the typical spectroscopic limit of i = 24.5; and sec- 
ondly, to reach at least to M* for galaxies up to z — 3. 

Throughout this paper, all magnitudes are given in the AB 
system unless stated otherwise. Unless otherwise stated, we use 
a ACDM cosmology given by Hq = 70 km/s/Mpc, Qo - 0.3 and 
A =0.7. 



2. Observations and data reductions 

2.1. Observations 

The NIR observation s were taken using the WIRCam detec- 
tor dPuget etal]|2004l) at the Canada France Hawaii Telescope 
(CFHT), with the exception of D2 J which was observed with 
WFCAM on UKIRT. WIRCam consists of four 2048 x 2048 
HgCdTe arrays arranged in a 2 x 2 format, with gaps of 45" 
between adjacent chips. The detector pixel scale is 0.3"/pixel 
resulting in a field-of-view of 21' x 21'. Observations were 
taken in a series of runs from 2005-2007 and were made in 
co-ordination with the COSMOS consortium. In this work we 



http://cadcwww.dao.nrc.ca/cfht/WIRDST0002.html 
http://terapix.iap.fr/article.php?id_article=832) 
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combine t he NIR imaging with th e CFHTLS T0006 imaging de- 
scribed by iGoranova et alJ (120091) . We note that in the D2 field, 
we include the COSMOS d ata that has already been described 
by iMcCracken et al.1 (1201 Ol) although we restrict ourselves only 
to the central 1 deg 2 and use this data in conjunction with the 
CFHTLS optical ugriz data (rather than the COSMOS Subaru 
observations) in order to be consistent between each of the four 
fields. 

Observation are summarised in Table Q] The mean exposure 
time listed is approximately the exposure time per pixel and is 
simply the individual exposure time multiplied by the number 
of exposures and divided by the number of subfields across the 
field. Our observing strategy was chosen in order to ensure uni- 
form coverage reaching a target depth of 24AB. Given the lim- 
ited observing time allocated to the program, this precluded full 
coverage of each of the 1° x 1° CFHTLS deep fields. Instead a 
subsection of each of the Dl, D3 and D4 fields was chosen (how- 
ever D2-COSMOS was covered in its entirety by the COSMOS 
consortium). Five WIRCam pointings were made in the Dl 
field and three in D3 and D4 fields. This gives a total area of 
0.6 deg 2 + 1.0 deg 2 +0.4 deg 2 + 0.4 deg 2 = 2.4 deg 2 (or 2.03 deg 2 
after masking). The subfields in Dl field were chosen to over- 
lap w ith the VIMOS-VLT Deep Survey (VVDS; iLe Fevre et al] 
2005) and the Spit zer Wide-area Infrared Ex tragalactic Legacy 
Survey (SWIRE; Lonsdale et al]|200l 12001 . whilst the entire 
field overlaps with a portion of the XM M Large Scale Structure 
Survey (XMM-LSS; IPierre et al.ll2 004). Two of the pointings in 
the D3 field were chosen to provide near-infrared imaging for a 
segment of the DEEP2 spectroscopic survey and the AEGIS X- 
ray survey. Finally, the three D4 pointings were chosen to coin- 
cide with a concurrent program to obtain spectroscopic redshifts 
of z ~ 4 galaxies using VLT VIMOS (PI: J. Bergeron, 077.A- 
0357). 

Observations were conducted using J, H and filters. 
Transmission plots of these filters are available from CFHT0. 
The integration times per exposure for all J, H and K s band ob- 
servations was 45s, 15s and 20s respectively. The observations 
were carried out in queue-scheduled mode with image quality 
constraints of 0.55" < FWHM < 0.65" and were "micro- 
dithered" (the telescope was displaced using sub-pixel offsets) 
using the standard WIRCam micro-dither pattern consisting of 
a 2x2 dither patter with offsets between consecutive dithers of 
0.5 pixels. Since WIRCam has a pixel scale of 0.3"/pixel, this 
micro-dithering is required in order to produce well-sampled 
images under our seeing constraint and also to allow matching 
with the CFHTLS pixel-scale of 0.186"/pixel. A further large- 
scale (« 1 - 2') dithering pattern was applied to the observations 
to avoid gaps in the coverage due to the gaps between adjacent 
CCDs. 

2.2. Data reduction 

2.2.1. WIRDS D1, D3 and D4 fields 

The WIRDS data set were reduced at CFHT and TERAPIX. 
Initial processing was made at CFHT using the I ' iwiQ prepro- 
cessing pipeline, which incorporates detrending and initial sky- 
subtraction processing. In order to ensure precise measurements 
of fluxes in background-dominated near-infrared observations 
like WIRDS, we use a two-step reduction process in which ob- 
jects were flagged in initial stacks created using a standard sky- 



3 cfht.hawaii.edu/Instruments/Filters/wircam.html 

4 cfht.hawaii.edu/Instruments/Imaging/WIRCam/IiwiVersionlDoc. 



subtraction procedure. These objects were identified and masked 
in a second-pass sky-subtraction. The overall reduction process 
was as follows: 

1 . Bias subtraction and flat-fielding (detrending) 

2. Initial sky-subtraction using object mask generation from the 
individual images 

3. Cross-talk correction 

4. Astrometric calibration and production of initial pre- 
processed stacks 

5. Second-pass sky-subtraction using object mask from the pre- 
processed stack and individual images 

6. Cross-talk correction 

7. Astrometric and photometric calibration 

8. Production of final processed image stacks 

These steps are now described in detail. The detrending 
stage involves the treatment of a number of instrumental im- 
prints: flagging saturated pixels, non-linearity correction, ref- 
erence pixel subtraction, dark-frame subtraction, dome flat- 
fielding, bad-pixel masking and guide window masking. The 
pre-processing of the images consists of bias-subtraction and 
flat-fielding using the bias and flat-field frames provided by the 
CFHT WIRCam queue observing team. A global bad pixel mask 
was generated using the flat to identify the dead pixels and the 
dark frames to identify hot pixels. For each image, we used the 
TERAPIX tool QualityFITS to produce weight-maps, object 
catalogues and overviews of individual image qualities (e.g. see- 
ing, depth). The productio n of the weight maps ma de use of the 
WeightWatcher software (Mar mo & Bertinll2008l) . 

After detrending, the initial sky-subtraction process was 
made. Median skies were produced for each individual image 
using up to 30 (object masked) images (equivalent to ~ lOmins 
total exposure time) in a given dithering pattern with on-sky sep- 
arations of < 30". The images were then sky-subtracted using 
these sky frames. Sky-flat images were constructed for each im- 
age using combinations of masked adjacent images taken over 
time intervals At and angular separations of A6 of Af < 20 mins 
and Ad < 10' respectively. In general individual sky-flat images 
consisted of « 8-9 multi-extension fits image files. To per- 
form the masking o f the images prior to co nstructing the sky- 
flats, SExtractor dBertin & Arnoutsl [l996) was used to iden- 
tify objects using a detection threshold of 1.5<x, a minimum de- 
tection area of four pixels, a background mesh size of 64 pixels 
and a background filter size of one pixel. A "check-image" of 
type "OBJECTS" was generated to identify objects in each im- 
age. With the images masked, the sky-flats were constructed and 
subtracted from the original images. 

After sky-subtraction, we cross-talk corrected our images. 
Cross-talk arises as a repetitive signal across the WIRCam am- 
plifiers on a given array and although cross-talk signals have 
been dealt with at the hardware level on WIRCam since 2008, 
many of our observations were taken prior to this date and so 
it must be dealt with in software. The cross-talk appears as 
"ghosts" of bright objects in adjacent amplifiers in a given ar- 
ray in the WIRCam detector. Depending on the form, the effect 
can be seen in all 32 amplifiers in a given array (where each am- 
plifier consists of 64 x 2048 pixels). We attempt to correct for 
the effect using a form of the "Medamp" procedure, where ob- 
ject masks were first created for each image using SExtractor. 
For each array, we then take the median of the 32 amplifiers 
as a way of identifying those features that are common to all 
32 amplifiers (i.e. the cross-talk signals). These signals are then 
(subtracted from the original image. 
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Table 1. Overview of the WIRDS/COSMOS observations. 



Field 


Subnelds 


R.A. 


Dec 


Area (masked) 


Band 


Ind. exp. 


N im 


Mean exp. time 






(J2000) 


( deg 2 ) 




(s) 




(s) 


CFHTLSD1 


5 


02:26:00 


-04:30:00 


0.6 (0.49) 


J 


45 


1,460 


13,140 










0.6 (0.49) 


H 


15 


4,916 


14,748 










0.6 (0.49) 


K s 


20 


4,235 


16,940 


CFHTLS D2/COSMOS 


9 


10:00:29 


+02:12:21 


1.0 (0.80) 


J 


5/10 


15,604/5,447 


13,710 










l.O (0.80) 


H 


15 


14,000 


13,125 










l.O (0.80) 


K s 


20 


7,272 


9,090 


CFHTLS D3/EGS 


3 


14:17:54 


+52:30:31 


0.4 (0.39) 


J 


45 


1,079 


16,185 










0.4 (0.39) 


H 


15 


3,390 


16,950 










0.4 (0.39) 


K s 


20 


2,622 


17,480 


CFHTLS D4/Q2215-1744 


3 


22:15:31 


-17:44:05.4 


0.4 (0.36) 


J 


45 


1,258 


18,870 










0.4 (0.36) 


H 


15 


3,166 


15,830 










0.4 (0.36) 


K s 


20 


2,362 


15,746 



Astrometric calibration was made using Scamp dBertinl2006l) 
with 2MASS as the astrometric reference catalogue. We ob- 
tained astrometric solutions with internal accuracies of « 0.1". 
Initial median-c ombined stacks were constructed using Swarp 
(iBertin et al.l20 02) at the instrument pixel scale. From the initial 
stacks we constructed a full object mask incorporating the fainter 
objects not detectable in individual images using the same time 
techniques for the individual object masks. Following this, the 
sky flat creation was repeated on the original detrended images, 
using both the individual image masks and the stack-mask to 
mask objects from the images before making sky-flats. The new 
sky flats were then subtracted from the detrended images. Cross- 
talk correction was again performed on the images. Astrometric 
calibration was performed using Scamp and the photometric cali- 
bration was then performed by matching to 2MASS photometry. 
All images were scaled to have a zero-point of 30.00AB. 

Before the final stack, we performed quality assessment of 
the individual images. First each image was inspected and those 
with severe defects were removed. We measured the seeing for 
each image based on the flux radius using PSFEx (this is defined 
as the radius which includes 50% of the flux with respect to the 
total flux) . All the final stacked images were combined after 
discarding the 5% worst-seeing in each field and filter combi- 
nation. The distributions of the flux_radius measurements for all 
images for each field/filter combination are shown in Fig. Q] In 
each case, the dash-dotted, dashed and solid lines line shows the 
flux_radius distribution across all JHK S respectively (as a frac- 
tion of the total number of images). Distributions are sharply 
peaked with a small tail of poor-seeing images (which are re- 
moved from our final sample by our 5% cut). The 5% cuts for 
each field/filter combination and the measured seeing for each 
final stacked image are given in Table [2] 



Table 2. Flux-radius properties of the WIRDS data. The seeing 
cuts give the upper limit above which individual exposures were 
rejected prior to stacking the images using Swarp. 



Field 




Seeing cut 




Final fluxjadius 




J 


H 


K s 


J 


H K s 


Dl 


0.90" 


0.87" 


0.91" 


0.68" 


0.62" 0.67" 


D2 




1.1" 


1.1" 


0.91" 


0.69" 0.68" 


D3 


0.82" 


0.81" 


0.89" 


0.64" 


0.58" 0.60" 


D4 


0.77" 


0.72" 


0.70" 


0.59" 


0.57" 0.58" 





0.25 




0.20 




0.15 




0.10 


c 




o 


0.05 


o 
2 


0.00 


LL. 


0.25 




0.20 




0.15 




0.10 




0.05 




0.00 



- D1 




J 


- D2 






H : 




/ 
ii 
ij 
tl 
i j 


\ 

vV* 
' \\ 

\ \ 


K — : 




1 1 








ij . 

1 1 1 «f if 




"^VTT-i i 




- D3 






- D4 

\. 




\ 






i 
f 


\ .-"•>. 






i j 
1 1 

'/ 

J ' 


s\ \ 
! V * 




i \ 


0.4 


0.6 0.8 


1.0 1.2 


0.4 0.6 0.8 1.0 1.2 



Seeing (") 

Fig. 1. Flux_radius distribution of the individual exposures in 
each of the four fields. The dashed-dotted, dashed and solid lines 
show JHK s -bands respectively (note that the COSMOS /-band 
data is not shown here, as we do not have access to the individual 
images used in the stack). 



Before stacking, images were resampled using a Lanczos-2 
4-tap filter with a 128 pixels mesh for background subtraction. 



Stacks were created in each band by sigma-combining appropri- 
ately weighted pixels using a modified version of Swarp kindly 
supplied by S. Foucaud; the pixel rejection limit was set to 3cr. 
Images were rescaled to match the pixel scale of 0.1 86" /pixel 
and tangent points of the CFHTLS deep stacks, making it easy 
to extract matched catalogues in double image mode. We note 
that at this point we performed a qualitative check on the as- 
trometry of the images by constructing combined colour i mages 
of the final stacks in each fields using Stiff package (Bertin 
l2010h . 



2.2.2. D2/COSMOS data reductions 

The COSMOS H band data was reduced identically as the other 
fields described above. However, the WIRCam K s and WFCAM 
/ band data were processed using a different pipeline which is 
described in lMcCracken et al.ld2010b The W FCAM J band data 
is partially described in lCapak et alJ d201 ll) . 
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WIRDS D1 , J WIRDS D1,H WIRDS D1, Ks 




150.6 150.4 150.2 150.0 149.8 150.6 150.4 150.2 150.0 149.8 150.6 150.4 150.2 150.0 149.8 

R.A, (degrees) R.A. (degrees! R.A, (degrees) 



Fig. 2. Weight maps for the Dl and D2 fields for each band. Each image covers 1° x 1°; black areas are regions without data. The 
grey-scale gives the 3<x depth for 5 pixel sources in 2" apertures. Note that crosses of slightly shallower depth are present in each 
pointing due to the gaps between detectors (broad central crosses in each pointing) and also due to the use of single guide stars for 
many observations (narrow crosses present in individual quadrants in some pointings). 



2.3. Completeness and coverage 

Each panel in Fig.[2](Dl and D2) and Fig.[3](D3 and D4) shows 
the weight maps for each field and covers 1 ° x 1 . The grey-scale 
represents the depth measured as the 3<x limit, based on a mini- 
mum object area of 5 pixels. This illustrates the excellent unifor- 
mity of the observations, especially in the final H and K s band 
stacks. We note that some depth variability between pointings 
can be seen in the /-band weight-maps in the Dl and D2 fields, 
which is due to loss of observing time through poor weather. 
Note that black regions in the thumbnails show areas lacking in 
any coverage. The COSMOS/D2 field is the only field with com- 
plete coverage of the entire CFHTLS Deep 1° x 1°. The remain- 
ing three fields cover (prior to any masking) of 0.6 deg 2 (Dl), 
0.4 deg 2 (D3) and 0.4 deg 2 (D4). 

We estimate the completeness of the final images by adding 
and detecting stellar sources in the stacks. Each image is split 
into individual pointings (5 for Dl, 9 for D2 and 3 for D3 and 
D4) of ~ 20'x20'. The 50% completeness limit for point-like ob- 
jects is then estimated by placing simulated stellar-like objects at 
a random position for a given pointing with the appropriate PSF. 
We then use SExtractor to detect the simulated objects and 
measure the percentage of such objects successfully recovered 
based on 500 objects for a given magnitude interval. 

Fig. [4] shows our results for all 12 final images. Detection 
success rates are plotted for each individual pointing separately; 
it can be seen that our fields have similar completeness limits. In 
each case, we give the mean 50% limit across all pointings for 
a given image (straight dashed lines in each plot) along with the 
standard deviation between those pointings. The 50% complete- 
ness limits in the Dl, D3 and D4 fields are consistently larger 



than 24.5 AB in all three bands, whilst the D2 data reaches depths 
of « 24 AB in the H and K s bands and 23.4AB in the /-band im- 
age. 

Consistent with the weight-maps, this figure shows that 
largest variation between pointings is for the Dl, D2 and D3 
/- band images, however these are still consistent to <x ~ 0.2 
(as opposed to cr < 0.1 for the H and K s images). 

2.4. Object extraction and photometry 

We have performed object extraction on the final NIR images in 
single and dual-image detection mode using SExtractor with 
a three pixel minimum detection area and a detection threshold 
of 1 .5cr. A background mesh size of 128 pixels was used for the 
background sub traction. Flu xes and magnitudes were calculated 
using the Kron (lKronlll980h -like total magnitudes (MAGJVUTO) 
and 2" (10.75 pixels) and 4.5" (20 pixels) diameter aperture 
magnitudes. Note that for the MAG_AUT0 estimate, the minimum 
allowable aperture was set to 2.5". As seeing is relatively homo- 
geneous across all the fields and bands and we do not perform 
any PSF homogenisation before catalogue extraction. 

For dual-mode extraction, we used gri "x 2 " K s band im- 
ages as detection images, producing two sets of catalogues. In 
both cases we ran SExtractor in dual-mode with the detection 
image and all eight available images (i.e. ugrizJHKs) in each 
field from the CFHTLS T0006 Deep and WIRDS. The gri X 2 
was constructed from the CFHTLS T0006 Deep gri images with 
Swarp in the "CHI2" combination mod^. This produces an op- 
timal combination of the input images where each pixel repre- 



5 http://terapix.iap.fr/cplt/T0006-doc.pdf 
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Fig. 3. As in Fig.|2] but for the D3 and D4 fields. 



sents the probability of it being drawn fr om the sky di stribution, 
based on a reduced^ 2 of the background! Szalay et al. 1999;). We 
compare the final photometry in our catalogues using 2MASS 
(which original individual exposures were calibrated, in Fig. [5] 
Small offsets are observed which for the most part are < 0.04. 

It is well known that SExtractor underestimates photomet- 
ric errors f or faint sources in combined images with correlated 
noise (e.g. lLabbe etal J 120031; iLeauthaud etaLll2007l) . We esti- 
mate the magnitude of this effect by measuring the variance in 
blank apertures in random apertures. First a 4000 x 4000 subim- 
age is extracted and aperture photometry is made with 2000 x 2" 
apertures placed on blank regions of sky (note objects in the im- 
age are masked during this procedure). The standard deviation of 
the flux over all apertures is calculated. Next, we measure aper- 
ture photometry on all the SExtracted objects in the 4000 x 4000 
pixel segment using 2" apertures. The correction factor is then 
taken as the ratio between the standard deviation of the blank 
sky apertures and the mean of the SExtractor errors for all 
the objects in that region. We repeat this process for a total of 
five 4000 x 4000 regions around each image and take the final 
correction to be the mean ratio across the image. This results 
in scaling of the magnitude errors of: f u = 2.10, f g = 2.66, 
f r = 3.06, ft = 3.08, f z = 2.76, fj = 2.43, f H = 2.61 and 
fx s = 2.49. The final errors on the magnitudes are given by the 
SExtractor magnitude errors multiplied directly by f x (note 
that these corrections have not been applied to the publicly re- 
leased catalogues). For the final (publicly released) catalogues, 
we do not apply any magnitude cut. 

We e stimate galactic dus t-corrections for individual objects 
using the lSchlegelet"ail (fl998) dust maps. Values of A were esti- 
mated for each filter using the analytic formula of ICardelli et al.l 
(Il989h . giving R(u) = 4.705, R(g) = 3.619, R(r) = 2.679, 
R(i) = 1.987, R(z) = 1.522, R(J) = 0.870, R(H) = 0.571 and 
R(K S ) = 0.366. 
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In order to exclude regions where the photometry may be 
unreliable, for example due to bright stars or at the field edges, 
or unavailable, i.e. outside the NIR coverage, we apply masks 
to our final catalogues in order to flag objects. In each field, 
these masks are a combination of the CFHTLS optical masks and 
masks based on the WIRDS coverage. Once these masks have 
been applied to the catalogues, the total usable area is 0.49 deg 2 
(Dl), 0.8 deg 2 (D2), 0.4 deg 2 (D3) and 0.4 deg 2 (D4). 

2.5. Photometric redshifts and galaxy properties 

With our unique combination of optical and NIR data (with 
which we can detect the 4000A break to z ~ 4) it is pos- 
sible to estimate photometric redshifts and stellar masses re- 
liably over a large redshift range. W e use the "Le Phare'0 
(lArnouts et al.ll2002t lllbert et~al1 l2006) package to measure pho- 
tometric redshifts and galaxy properti es with a y 2 temp late- 
fitting meth od, based on the method of Ilb ert et a U (120061) and 
lllbert et all d2009l). As in lllbert et alJ d2009l) . we use a combi- 
nation of the Polletta et al J (I2007h templates with the additional 
templates generated by lllbert et al.l (|2009) to account for star- 
bursts with ages of 0.3 to 3 Gyr. Following lllbert et alJ {2010), 
we compute the final photometric redshifts using the median of 
the probability distribution function (PDFz). 

We use an extensive collection of spectroscopic redshifts 
to calib rate our photometric redshifts. In Dl, we use VVDS 
"Deep " dLe Fevre et alj |2005) and "Ultra Deep" (Cuccia ti et al.l 
2012; Le Fevre et al. In Prep) spectroscopic samples. The 
publicly-available VVDS Deep sample contains 8,981 objects 
over 0.5 deg 2 in Dl. It is a purely magnitude-limited sample 
with i < 24 and spans a redshift range of < z < 5. The 
Ultra Deep sample comprises ~ 1500 spectra over an area of 



cfht . hawaii . edu/~ arnouts/LEPH ARE/cfht_lephare 
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Fig. 4. Point-source object completeness for the J (left-hand panels), H (middle panels) and K s (right-hand panels) observations in 
the WIRDS fields. For each field/filter combination, the completeness is measured for each of the individual pointing regions (i.e. 
pl-p5 in the Dl, pl-p9 in the D2, pl-p3 in the D3 and pl-p3 in the D4), illustrating the uniformity of the fields. The dashed lines in 
each case indicate the mean 50% completeness for the given field/filter combination. 



« 0.14 deg 2 and covers a magnitude range of 22.5 < i < 24.75. 
Both VVDS spectroscopic catalogues contain redshift quality 
flags which range from 1 to 4 with 1 being most unreliable and 
4 being most reliable; in addition a flag 9 is given to objects 
identified based on a single emission line. In this analysis we 
only use quality flags 3 and 4 from the VVDS catalogue (where 
these are the objects with the most secure spectroscopic classifi- 
cations). Using the VVDS Ultra Deep data we find an outlier rate 
of T] — 3.7% and o"Az/(i +z ) = 0.025, with a median magnitude of 
i* ,. = 24.0, whilst using the Deep data we find an outlier rate 

median ' 1 

of 77 = 1.9% and cr A ,/ (1+z) = 0.030, with a median magnitude of 
i* , — 23.7. The median offset between the photometric and 

median 1 

spectroscopic redshifts is Az/(1 + z) = 0.017. 

In the D2/COSM OS field we use the zCOSMOS "10K" 
spectroscopic sample dLillv et al.ll2009h . We find 3004 objects 
predominantly in the magnitude range 17.5 < i < 22.5 and 
over a redshift range up to z < 1.4 present in our photomet- 
ric catalogue. From this data (and as before, using only objects 
with spectroscopic flags of 3 or 4) we estimate an outlier rate of 
r] = 1.4% and CTAzf(i+z) - 0.023, based on a sample with me- 
dian magnitude of ;' median = 21.6 (spectroscopic observations for 
the VVDS and zCOZMOS in Dl and D2 fields were made us- 
ing VIMOS on the VLT UT3). The median offset between the 
photometric and spectroscopic redshifts is Az/(1 + z) = 0.020. 

In D3, we use t he DEEP2 DR3 redshift catalogue 
dDavis et al.ll2003ll2007l) . based on observations using the Deep 
Imaging Multi-Object Spectrograph (DEIMOS) on Keck II. The 



catalogue contains 47,700 unique objects, of which 2,977 have a 
'zquality' flag of > 3 (i.e. are deemed to be reliable redshifts) and 
are present in our photometric catalogue. This sample predom- 
inantly covers a magnitude range of 1 8 < i < 24 with the bulk 
being below a redshift of z < 1.6. Using the DEEP2 data, we 
estimate an outlier rate of 77 = 3.0% and 0"Aj/(i+ Z ) = 0.027, based 
on a sample with median magnitude of * median = 22.5 for the 
WIRDS D3 photometric catalogue. The median offset between 
the photometric and spectroscopic redshifts is Az/(l+z) = 0.017. 

In the D4 field, we use spectra obtained using the AAOmega 
instrument on the Anglo-Australian Telescope (AAT) as part 
of a program to provide optical spectrosco py of X-ray point- 
sources in the CFHTLS (IStalin et al.l l2Q10h . The observations 
provide redshifts for 1,800 objects in the D4 field, of which 
168 are QSOs, 66 ar e stars and 1 , 335 are galaxies, all at mag- 
nitudes of i < 22.5 dBielbv etalJl2010h . In total, 1,090 of the 
galaxies overlap with our photometric data, most of which are at 
z < 0.8. Based on these, we find an outlier rate of 77 = 2.1% and 
o"Az/(i+z) = 0.021, with a median magnitude of * median = 20.0. 
The median offset between the photometric and spectroscopic 
redshifts is Az/(1 +z) = 0.013. 

In Fig. [6] photometric redshifts are plotted against spectro- 
scopic redshifts. The dashed line shows the outlier constraint of 
Az/(1 + z) = 0.15. In addition, we plot the redshift accuracy, 
o"Az/(i+z) as a function of /-band magnitude for all four fields 
combined in Fig. [7] These figures show that our photometric red- 
shifts are highly accurate up to at least z = 1.5 and i - 25, with 
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Fig. 5. Photometric offsets for stars in WIRDS compared to 
2MASS for each field and filter combination (note that these off- 
sets have not been applied to the publicly released catalogues). 



these limits imposed more by the available spectroscopic data 
than any reduction in photometric redshift accuracy. 



2.6. Estimating stellar masses 



Following lllbert et al.l (|2010) we use stellar population synthe- 
sis (SPS) models to estimate stellar masses effectively con- 




0.5 1.0 1.5 2.0 
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Fig. 6. Photometric redshifts versus spectroscopic redshifts for 
all four fields. Top left shows the WIRDS Dl data with VVDS 
Deep and Ultra-Deep spectroscopic redshifts (3,192 galaxies). 
The top right panel shows the WIRDS D2 data compared to 
zCOSMOS 10K redshifts (3004 galaxies). Bottom left shows the 
WIRDS D3 data with DEEP2 spectroscopic data (2977 galaxies) 
and bottom right shows the WIRDS D4 data with the AAOmega 
spectroscopic redshifts (700 galaxies). 



verti ng galaxy lum inosity to a stellar mass dBell et al.l 120031: 
iFontana et al.ll2004l) . 

The SED templates were generated with the SPS package 
developed by iBnizual & Charlotl (120031) . Although these SPS 
models do not include the contribution of stars in the ther- 
mally pulsing asymptotic giant branch (TP-AGB) phase, we note 
that recent studies have shown the contribution of this phase to 
have been overestimated (by factors of ~ 3) by models that do 
include it (iKriek et al.l l2010: Zib etti et al.1 12012). Significantly, 
IZibetti et al.l (120121) find that the NIR SEDs of post-starburst 
galaxies (a population that would be expected to exhibit the 
strongest influence from the TP-AGB phase) are consistent with 
the SEDs predicted by the original iBruzual & Charlotl (120031) 
models. 

We assumed a universal Chabrii^ (120031) IMF and an ex- 
ponentially declining star formation history of the form SFR 
oc e~^ r , with t in the range 0. 1 Gyr - 30 Gyr. The SEDs were gen- 
erated for a grid of 5 1 ages in the range 0. 1 Gyr - 1 4.5 Gyr. Dust 
extinc tion was applied to the templates using the ICalzetti et al.l 
(120001) law for < E(B - V) < 0.5. The templates were calcu- 
lated based on two different values for the metallicity, 0.008Z o 
and O.O2Z . We impose a prior that significant extinction is only 
allowed for galaxies with high SFR, such that E(B - V) < 0.15 
for galaxy ages of t > At. 

3. Galaxy populations and number counts 

3.1. BzK photometric selections 

One of the key scientific objects of the WIRDS project is to 
probe the galaxy population at z > 1.4, in particular via the 
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Fig. 7. The redshift accuracy, <ta z /(i +C ), as a function of /-band 
magnitude for all four fields (connected diamond points). The 
spectroscopic data used in this analysis are described in the text, 
whilst the numbers of galaxies used for the comparison of spec- 
troscopic versus photometric redshift at each magnitude is given 
by the grey histogram. 



"BzK" selection iDaddi etal] d2004l) . which is a means to isolate 
stars and galax ies of different typ es and redshifts. In the their 
original paper, IDaddi et"ail |2004) used BzK colours to iden- 
tify selection criteria to isolate galaxies at 1.4 < z < 2.5 and 
to separate galaxies within this range into likely star-forming 
and passive galaxy samples. This has since been used by many 
authors to identify such galaxies and inve stigat e their prop- 
erties. At lower redshifts, lLane et al . (200% and lHartlev et al. 
d2010l) is olate features of re presenting passively evolving galax- 
ies, with lLane et al.l d2007l) identifying "branch" galaxies and 
lHartlev et alJ d20 101) identifying star-for mation limits in the BzK 
colour space. However, as pointed out by Rigucc ini et al.l d201 ll) . 
colour selections can be potentially biased, wi th the effect of 
dust o ften significant. Using Spitzer 24 fim data, Rigu ccini et alJ 
d201lh measured the success rate of the BzK selection in identi- 
fying 1.4 < z < 2.5 star-forming galaxies and fo und it to be ~ 
90% c omplete to their magnitude limit. Similarly, Grazian et al. 
(2007) presented an analysis of the BzK selection, evaluating 
the limitations of the selection based on multi-wavelength data. 
They conlcude that the BzK selection is highly efficient at iden- 
tifying galaxies in the redshift range 1.4 < z < 2.5, but that for 
galaxies faint in the K band and red in the z - K colour, it is dif- 
ficult to distinguish between star-forming and evolved galaxies, 
resulting in an underestimation of the passively evolving popu- 
lation. 

In this Section we will present transformations from the gzK 
photometry of CFHTLS/WIRDS to the BzK colours of Daddi et 
al. and compare the selection results to that obtained based on 
our 8-band photometric redshift catalogue. 

In Figure [8] the solid blue, red and black curves show the B, 
z and K transmission profiles of the filters used by Dadd i et alJ 
(2004) and the dashed green, red and black curves show the 
MegaPrime g, z and WIRCam K s filters profiles. Although the 
z and K filters are similar, the MegaCam g filter is offset from 
the Subaru B filter. 
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Fig. 8. The gzK and BzK filter sets. The dashed lines show g 
(green line) and z (orange line) filters and the K s filter (black 
line). The solid lines show the Daddi et al, filters: Bsesseii (blue 
line), Zcunn (orange-line) and the ISAAC K s filter (black line). 
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Fig. 9. Comparison between COSMOS BzK from 
McCr acken et al.l (l2010l) and transformed WIRDS gzK colours 
for objects common between the two datasets. Orange and cyan 
points are stars galaxies respectively. 



We note that BzK colours for objects the D2/C O SMOS 
field has already been published McCracke n et alJ ([2010). 
iMcCracken et alJ ma ke a transformation from their B, z and 
K filters to those of Daddi et al. (2004) colours using stellar 
spectra. However, the McCracke n et al.l (120 1 Oh filters are over- 
all closer to the original IDaddi et all (120041) filters. In particu- 
lar, both works use the same B filter. Given the differences be- 
tween the B and g filters, the simple stellar spectra m ethod can- 
not be used to transform between our system and the lDaddi et al.l 
system. We therefore use the available COSMOS BzK colours 
for galaxies and stars in the D2 field to compute an empirical 
transformation and to account for differences in galaxy spectral 
slopes. We first re-derive the colour transformation using a linear 
form for (g - z)wirds to (B - z)Daddi based purely on the stellar 
component and for the (z - A^wirds to (z - A")Daddi using both 
the stellar and galaxy components. We then correct the stellar 
transformation of the (g - z)wirds to match the {B - z)Daddi as a 
function of the (u - g)wiRDS colours for the galactic component. 
The resulting offsets in the colours compared to the COSMOS 
data for both stars (orange points) and galaxies (blue points) 
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are shown in Fig. [9] The transformations used to convert the 
CFHTLS/WIRDS colours to the BzK colours of iDaddi et all are 
then given by: 



(B - zbaddi = 0.410 + 1.2 17(£ - z) wirds + C ug 
(z - K) DRddi = 0.033 + 0.987(z - /Owirds 



(1) 
(2) 



where C ug is the correction based on the (u - g) colour and is 
equal to: 



(3) 



C ug = -0.318 + 0.291O-g) 
or equal to zero if the above is greater than zero. 

3.2. Galaxy number counts 

We now present the galaxy number counts 
our fields. Galaxy counts are the basic 
galaxy population and provide an important 
test of any survey. Furthermore, near-infrared selected 
galaxy counts provide some are able to discriminate be- 



for each of 
and statistic 
zeroth-order" 



In order to separate s tars and galaxies we use the BzK colour 
cut given bv lDaddi et ali 
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Fig. 10. J (AB) number counts from each of the WIRDS fields. 
The filled circles, stars, triangles and squares give the galaxy 
number counts for the D 1 , D2, D3 and D4 fields respectively. For 
comp arison the red open ups i de dow n diamonds show the counts 
of ICristobal-Hornil los et al] d2009l) . the red asterisks show the 
2005) data and the red x's show the counts of 
1. 



lovino et al 



Saracco et al "J (1200 1 



(z - K s ) > 0.3(B - z) - 0.5 



(4) 



Before applying this cut, we transform our colours to BzK 
as described above and apply the photometric offsets given in 
Fig. [5] The resulting galaxy number counts in each of the three 
NIR bands are given in Table[3] combined across the four fields. 
In addition, we also give the standard deviation between the four 
fields for the counts in each band. 

The number counts of galaxies in each of the four fields 
are shown in Fig. [10] (/-band), Fig. QT| (//-band) and Fig. [12] 
(Zf s -band). These are all measured directly from the individ- 
ual J, H or K s selected catalogue in each case (as opposed 
gri- or K s - selected catalogues). In each case Dl, D2, D3 and 
D4 counts are given by circles, stars, triangles and squares re- 
spectively. Both galaxy counts (filled black points) and stel- 
lar counts (open orange points) are shown for each field, 
with literature counts are plotted for comparison. For the J 
band co unts we plot count s from the VIRMOS Deep Imaging 
Surv ey (llovino et all 120051 diamond s), the ALHAMBRA sur- 
vey (ICristobal-Hornillos et alJ I2009L asterisks) and combined 
deep-NIR counts fro m the Chandra Deep and HDF South fields 
dSaracco et al. I l200ll x's). In all four fields our number counts 
are consistent with those from the literature to J a 23 beyond 
which our galaxy counts begin to be affected by incompleteness. 

In //-band, we show the WHDF counts of [Metcalfe et al.l 
(120061 diamond s) and the HST/NICMOS observations of 
lYan et all (1 19981 asterisks). Again, close agreement is seen be- 
tween the WIRDS number counts and the previously published 
data, with the completeness beginning to affect the WIRDS 
counts at H m 23.5. 

Finall y, the K s band counts are shown in comparison with the 
counts of McCrack en et al.l (1201 Oi diamonds) in the COSMOS 
field (of which our data in cludes a subs e t), the VIMOS Deep 
Imaging Survey counts of llovino et~ail (120051 asterisks) and 
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Fig. 11. H (AB) number counts from each of the WIRDS fields. 
The filled circles, stars, triangles and squares give the galaxy 
number counts for the Dl, D2, D3 and D4 fields respectively. 
For comparison the red open diamonds show the counts of 
iMetcalfe et al.l d2006l) and the red asterisks show the lYan et all 
(119981) data. 



the Chandra Deep a nd HDF South fields combined counts of 
ISaracco etall(l200ll x's). 

Our stellar number counts show significant field-to-field 
variation, representing the different stellar populations in each 
deep field, with D4 having the highest density of stars and the 
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Table 3. Differential number counts of galaxies in the J, H and K s bands. 
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(/0.5mag/deg 2 ) 




(/0.5mag/deg 2 ) 




(/0.5mag/deg 2 ) 


15.75 


5 


2.5 


1.2 


40 


19.7 


3.8 


49 


24.1 


2.1 


16.25 


74 


36.4 


1.7 


107 


52.7 


5.3 


119 


58.6 


6.0 


16.75 


149 


73.3 


7.9 


213 


104.8 


12.7 


252 


124.0 


12.5 


17.25 


296 


145.7 


14.6 


397 


195.4 


24.6 


478 


235.2 


29.9 


17.75 


491 


241.6 


27.8 


721 


354.8 


46.7 


848 


417.3 


43.4 


18.25 


849 


417.8 


46.4 


1169 


575.3 


72.7 


1445 


711.1 


89.4 


18.75 


1436 


706.7 


87.8 


2000 


984.3 


151.3 


2608 


1283.5 


166.3 


19.25 


2393 


1177.7 


145.0 


3347 


1647.1 


215.3 


4340 


2135.8 


252.8 


19.75 


3880 


1909.4 


224.0 


5282 


2599.4 


368.5 


6842 


3367.1 


411.1 


20.25 


6133 


3018.2 


354.0 


8057 


3965.1 


493.3 


10183 


5011.3 


589.3 


20.75 


8986 


4422.2 


466.6 


11846 


5829.7 


708.0 


14069 


6923.7 


796.8 


21.25 


12965 


6380.4 


613.2 


16407 


8074.3 


934.9 


19337 


9516.2 


1030.1 


21.75 


18634 


9170.3 


830.7 


22666 


11154.5 


1130.0 


26239 


12912.9 


1340.0 


22.25 


26433 


13008.4 


1260.6 


30417 


14969.0 


1437.5 


35093 


17270.2 


1771.0 


22.75 


36960 


18189.0 


1994.7 


38927 


19157.0 


1699.8 


44340 


21820.9 


2173.7 


23.25 


45692 


22486.2 


2456.2 


47804 


23525.6 


2133.9 


53912 


26531.5 


2419.8 


23.75 


48400 


23818.9 


1888.6 


50023 


24617.6 


2033.5 


58358 


28719.5 


1924.6 
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Fig. 12. The K s (AB) number counts from each of the WIRDS 
fields. The filled circles, stars, triangles and squares give the 
galaxy number counts for the Dl, D2, D3 and D4 fields respec- 
tively. For comparison th e red open diamonds show the counts of 
McCrack en et alj yOlO), the red asterisks show the lovino et al 



(2005) data and the red x's show the counts of ISaracco et al 

hooib . 



Dl and D3 having lower densities. This is consistent with the 
galactic latitudes of the fields. 

3.3. Selection of star-forming and passive galaxies 

In this Section, we investigate a range of galaxy properties as 
function of type, i.e. star-forming versus passive galaxy popu- 
lations. There are a number of ways to accomplish this, for in- 
stance by using observed colours, rest-frame colours or full tem- 
plate fitting. For the remainder of this paper, we will focus on 
using the derived rest-frame colours as a consistent method to 
separate star-forming and passive galaxy populations. For this, 



we follow llibert et al.1 ((2010) and perform a selection in the de- 
rived rest fram e colour Mnuv - Affl (after dust correction is 
performed). In llibert et all (|2010), three selection criteria are 
used to identify star-forming, passive and intermediary popu- 
lations. Here we simplify this to two populations, star-forming 
(incorporating t he star-forming and intermediary populations of 
Ilber t et al.l2010l) and passive. Thus the selection criteria take the 
following form: 



M NUV - M r > 3.5 (Passive) 
M N uv - M r < 3.5 (Star - forming) 



(5) 
(6) 



The separation of star-forming versus passive populations in 
this way is proven to be a good indication of star-f ormation ac- 
tivity (e.g. iMartin et all 120071; lArnouts et"ai],2007l) and its suc- 
cess has been hi ghlighted by the co mparison with morphological 
classifications in llibert et al.l (l2010l) . 

3.4. Galaxies atz< 1.4 

We first look at the properties of the z < 1 .4 galaxy population 
in the WIRDS data, making use of the full grix 2 selected photo- 
metric catalogue. Fig.[T3lshows the BzK diagram for all stars and 
galaxies in WIRDS with photometric redshifts of z < 1 .4 and a 
magnitude limit of K s < 22. Orange points show the stellar pop- 
ulation (identified using the BzK cut), red the passive galaxies 
and blue points show the star-forming galaxy population (identi- 
fied as described above using the rest-frame Mnuv ~ M r colours 
of the gal axies). The dashed lines show the BzK colour cuts pre- 
sented by Dad di et all ([2004), which are given by equation 3| for 
the stellar-galaxy separation, whilst the selection of star-forming 
galaxies at z > 1 .4 is given by: 

(z - K s ) >(B-z)- 0.2 (7) 
and the boundary for passive galaxies at z > 1 .4 is given by: 
(z - K s ) <(B-z)- 0.2 n (z - K s ) > 2.5 (8) 



7 Here M NUV is based on the GALEX rest-frame NUV filter and M T 
is based on the SDSS r band filter. 
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Fig. 13. The BzK diagram for galaxies at z < 1 .4. The galaxy 
population is split into two groups: red points for passive galax- 
ies, and blue points for star-forming galaxies (separated using 
rest-frame Mnuv - M r colours). Only objects with K s < 22 and 
a photometric redshift of z < 1-4 are shown. T he black dashed 
lines show the cuts used by iDaddi et alj (|2004) to select z > 1 .4 
galaxies and stars, whilst the solid black line shows the model 
track for a passively evolving elliptical galaxy from z - 1 .5 to 
z — 0. The filled black circles on this line are at intervals of 
Az = 0.25, beginning at z = 0.25 and finishing at z = 1.5. 



These selection criteria follow closely the boundaries of our 
BzK transformed colour selection and photometric redshift lim- 
its to a good degree, with few z < 1.4 galaxies lying outside 
the z > 1.4 criteria. Based on our photometric redshifts, we find 
that only 0.5% and 0.5% of all K < 24, z < 1.4 galaxies fall 
within the z > 1.4 star-forming and passive galaxy selections 
respect ively. The st e llar l ocus also appears to be well separated 
by the IDaddi etafl (120041) criteria, suggesting th at our gzK to 
BzK tr ansformation are correct. As discussed by lHartlev et al.1 
(2010), there is a separation of star-forming and passive galaxies 
in the BzK at z < 1 .4 as well as at z > 1 .4. The passive galaxies 
in fact form a locus approximatel y parallel to the stellar locus 
at B - z > 3, already noted bv lLane et all (120071) in UKIDSS 
data. We highlighted this in Fig. [13] with the evolutionary track 
(black line with filled circles at intervals of Az = 0.25) of a typ- 
ical passively evolving galaxy fro m z = 1 -5 to z = 0, calcula ted 
using the galaxy evolution code of Bruzual & Chariot (2003), a 
Salpeter IMF and an instantaneous starburst and a redshift of 
formation of zj — 5. 

3.5. Passive and star-forming galaxies atz> 1.4 

We now turn to galaxies at z > 1.4, plotting once again the 
BzK diagram for the WIRDS data in Fig. [14] All points from 
Fig.[13]are plotted with grey points (i.e. stars and z < 1.4 galax- 
ies). Overplotted on this z < 1.4 population are galaxies with 
z > 1 .4 with blue points showing star-forming galaxies and red 
points showing passive g alaxies. Again the dashed lines show 
the colour constraints of IDaddi et all (120041) as given in equa- 
tions |4] [7] and [8] The sBzK constraint corresponds closely to 
galaxies with z p h t > 1 .4 and with high star-formation rate. At 



Fig. 14. The BzK diagram for z > 1.4 galaxies. The blue and 
red points show star-forming and passive galaxies at z > 1.4 
(based on the Mnuv - M r rest-frame colours and photometric 
redshifts). The grey points show the stellar and z < 1.4 galaxy 
populations (as plo tted in Fig. T3). Th e dashed lines show the 
selection criteria o f IDaddi et al.U2004 . 



K < 24, 23% of objects are scattered into the z < 1.4 region, 
however the bulk (76%) are at z > 1.4. The final 1% of star- 
forming objects (based on the rest-frame colour selection) are 
found to fall within the z > 1 .4 pass ive galaxy selection. W e find 
a somewhat lower success rate than Ri guccini et al.1 (1201 ll) . who 
found a success rate of 90%. This is probably a consequence 
of the increased scatter in our photometric redshift estimates in 
comparison with their work, which includes many more photo- 
metric bands in addition to, crucially, Spitzer 24 /im measure- 
ments in their photometric redshifts. However, both results ulti- 
mately show the effectiveness of the BzK selection in identifying 
star-forming galaxies at z > 1 .4. 

Concerning the passive z > 1.4 population, passive objects 
selected from the photometric redshift catalogue significantly 
overlap with the z < 1 .4 population. Many passive z > 1 .4 galax- 
ies are missed by the pBzK cut when applied to our data, falling 
in either the z < 1.4 region (41%) or the star-forming z > 1.4 
region (7%). In addition, the pBzK selection includes a signifi- 
cant fraction of z < 1.4 objects. We note that the pBzK selection 
relies on extremely deep B (or in this case g) band data and given 
the g band depths of « 26, these numbers for the pBzK selection 
are likely affected by incompleteness and increased photometric 
errors compared to deeper datasets. Despit e this, we note tha t 
similar results have been published, with iGrazian et al.1 (120071) 
for example reporting a completeness of only 34% for the pBzK 
selection in terms of retrieving the passive 1.4 < z < 2.5 galaxy 
population. 

Fig. [15] shows the photometric redshift distribution for at 
K s < 24 for pBzK and sBzK selected galaxy populations. The 
top panel shows the entire K s < 24 galaxy redshift distribution 
(solid black line). This extends well out to z ~ 3.5 - 4. 

The middle panel shows the sBzK galaxy redshift distribu- 
tion (solid blue histogram). We also show the z > 1.4 star- 
forming population as selected using the photometric redshift 
catalogue in which star-forming galaxies are identified via the 
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Fig. 15. Top: Photometric redshift distributions of the entire 
galaxy population at K s < 24. Middle: Redshift distribution of 
galaxies selected to be star-forming at z > 1.4 via the sBzK se- 
lection (solid blue histogram) and using photometric redshifts 
and template fits (dashed blue histogram). The dotted line shows 
the entire K s population. Bottom: Redshift distribution of galax- 
ies selected to be passive at z > 1 .4 via the pBzK selection (solid 
red histogram) and again using photometric redshifts and tem- 
plate fits (dashed red histogram). 



rest-frame colours (dashed blue histogram). Comparing the two, 
the photometric redshifts suggest that the sBzK selection misses 
a small fraction of the z * 1.4 galaxy population, whilst includ- 
ing a small number of galaxies in the range 1 < z < 1 .4. 

In the bottom panel of Fig. Q3J we show the pBzK redshift 
distribution (solid red histogram). We also show the z > 1.4 pas- 
sive population identified in the template fitting using Le Phare 
(dashed red histogram). In agreement with the BzK distribution 
we observe in Fig. [14] we see a significant contribution to the 
pBzK selection from galaxies identified with photometric red- 
shifts of 1 < z < 1.4. 

In form, the BzK redshift distri b utions are similar t o 
those measured by iReddv et alJ (l2005 | ). [G razian et al] (120071) . 
lHartlev et ai](l2008l) and lMcCracken et alJd2010l) , with the sBzK 
selection showing a relatively sharp lower-redshift cut-off at 
Z = 1 .4 and the pBzK selection extending to lower redshifts of 
z ~ 1 . Additionally, the sBzK selection extends to relatively large 
redshifts of z > 3. In contrast, few pBzK galaxies are observed 
with photometric redshifts of z > 2. 

We show number counts for sBzK samples in Fig. IT6l Also 



shown are counts of iKong et al.l (120061 triangles ! iBlanc et al 
(2008, diamonds) from the MUSY C Survey. lHartlev et al 
(2008, asterisks) from UKIDSS and IMcCracken et al.1 d2010l 
crosses) from COSMOS. In general we see good agreement be- 
tween our counts of sB zK galaxies a n d the literature counts, 
alt hough the counts of lHartlev et al.l (120081) and in particu- 
lar IBlanc et all (l2008h are marginally higher than the WIRDS 
counts. The sBzK number counts combined across all four 
WIRDS fields are given in Table g] 
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Fig. 16. Number counts of galaxies selected using the sBzK 
colour selection. Blue circles, stars, triangles and squares show 
the number counts of sBzK galaxies in the WIRDS Dl, D2, 
D3 and D4 fields respectively. Black points show compari- 
son coun ts from the literat ure, with the triangles showing the 
counts of lKong et all (12006), diamonds the counts of Blanc et al. 
(2008 1), asterisks those of [Hartley et all d2008l) and crosses those 
of IMc Cracken et all d2010h . 

Table 4. WIRDS sBzK number counts in the four WIRDS field 



K S (AB) 






(AB) 




(/0.5mag/deg 2 ) 


18.25 


3 


1.5 


18.75 


6 


3.0 


19.25 


22 


10.8 


19.75 


58 


28.5 


20.25 


278 


136.8 


20.75 


892 


439.0 


21.25 


2385 


1173.7 


21.75 


5169 


2543.8 


22.25 


9007 


4432.6 


22.75 


11155 


5489.7 


23.25 


6332 


3116.1 


23.75 


638 


314.0 



The WIRDS number counts for the pBzK samples are 
shown in Fig. [17] The literature counts are f r om th e same 
sources as given in Fig. [161 with iKong et al.l {2006) counts 
given by triangles, IBlanc et al.l (120081) counts denoted by dia- 
monds. |Hartlev et aLl(l2008l) given by asterisks and the counts of 
IMcCracken et alJ y010) given by crosses. Once again, the in- 
dividual WIRDS fields are consistent with each other and are 
in good agreem e nt wit h literature counts, except in the case 
of lHartlev et al.l ([2008), which are somewhat lower than the 
other pBzK counts. The WIRDS pBzK number counts combined 
across all four W IRDS fields are given in T able [5] 

As noted by McCracken et al. (2010), a turn-over in the 
pBzK number counts is seen at K s > 21.5, which indicate a 
lack of passive galaxies at z > 2. We note also that we see few 
passive galaxies at z > 1.8 identified by rest-frame colours in 
our photometric redshift sample (see Fig. [TBI . This is partially 
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Fig. 17. WIRDS number counts for pBzK galaxies in the D 1 (red 
circles), D2 (red stars), D3 (red triangles) and D4 fields (red 
squares). Literature counts are from the same sources and using 
the same symbols as Fig. [16] 



Table 5. WIRDS pBzK number counts in the four WIRDS field 



K s 




iVgal 


(AB) 




(/0.5mag/deg 2 ) 


19.25 


20 


9.8 


19.75 


91 


44.8 


20.25 


355 


174.7 


20.75 


699 


344.0 


21.25 


964 


474.4 


21.75 


814 


400.6 


22.25 


412 


202.8 


22.75 


111 


54.6 




Fig. 18. Estimated mass limits for the photometric catalogues in 
each of the four WIRDS fields. The greyscale contours show 
the distribution of the i < 25.5 galaxy population, normalised by 
area. Blue triangles with horizontal error bars show the estimated 
mass limits for consecutive bins in redshift corresponding to a 
magnitude cut of i = 25.5 for all galaxies. 

Table 6. Mass limits based on a maximum 30% of objects at 
magnitudes of i > 25.5 



Redshift 


Mass Limit 




(log 10 (M/M Q )) 


0.2 < z < 0.4 


7.75 


0.4 < z < 0.6 


8.26 


0.6 < z < 0.8 


8.31 


0.8 < z < 1.0 


8.67 


1.0 < z < 1.2 


9.98 


1.2 <z< 1.5 


9.65 


1.5<z<2.0 


10.61 



due to the limiting magnitude of our optical data. For a typical 
passive galaxy, we expect a z - K s ~ 2.5 and B - z ~ 4, which 
for K s w 21.5 would suggest optical magnitudes of g ~ 28 and 
z ~ 24. This is clearly a challenge for our g-band data, which 
have mean 50% completion limits of g = 26.2 / g = 26.7 for 
extended/point-like objects in each of the four fields. 



4. Galaxy mass functions 

We now present an analysis of the galaxy mass function based 
on the WIRDS photometric gn';^ 2 -selected catalogues. The stel- 
lar mass function traces the build-up of the stellar mass con- 
tent of galaxies, which provides a key observational constraint 
in modelling galaxy evolution and the physics that regulates it. 
With the deep NIR data and 2.4 deg 2 coverage over four fields, 
WIRDS provides an excellent dataset for such analysis, allow- 
ing the measurement of the mass function over a broad range of 
stellar mass, whilst the survey strategy significantly reduces the 
impact of cosmic variance on the results. 



4.1. Estimating the WIRDS mass limits 

Before calculating the stellar mass functions of our galaxy popu- 
lations, we first evaluate the mass limits of the galaxy population 
in the four WIRDS field s as a function of r edshift. This is done 
following the method of Ilbert et al. (2010), in which mass lim- 
its are calculated as the lowest mass at which less than 30% of 
galaxies are fainter than a chosen magnitude limit. We estimate 
this mass limit as a function of galaxy redshift and type (i.e. pas- 
sive, star-forming and both combined) using a magnitude limit 
of/ = 25.5. 

Fig. n~8l shows the resulting mass limits as a function of red- 
shift for each field. The grey-scale contours show the galaxy 
population distribution at i < 25.5 normalised by area. Mass 
limits for the galaxy population are shown by the blue triangles. 
The extent of each redshift bin is given by the horizontal error 
bars. Estimated limits across the four fields are broadly consis- 
tent based on the imposed i < 25.5 limit. In each of our five 
redshift bins we have taken the median mass limit of the four 
fields. These are given in Table|6] 
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Fig. 19. The mass function is shown as a function of redshift in each of the four WIRDS fields (blue circles, green diamonds, yellow 
triangles and red squares for the Dl, D2, D3 and D4 fields respectively). Measurements in each field were performed using different 
mass bins in order to be presented clearly in the plotting (i.e. there is no artificial offset between the plotted mass functions). Errors 
on the points are a combination of the Poissonian errors with the uncertainties from the photometric redshift estimates. The double 
Schechter fits to the data are shown by the long-dashed curves, whilst the short-dashed line that shows the 0.2 < z < 0.4 WIRDS 
mass function replicated in each panel for comparison. 



4.2. The stellar mass function from z = 0.2 to z = 2 

We calculate the mass functions using t he ALF (Algorithm for 
Luminosity Function; lllbert et all 120051) tool with a step-wise 
maximum likelihood (SWML; Efsta thiou et al.lfl988l) estima- 
tor. We note that although we only present the SWML results 
here, the mass function s have also be en c alculated using the 
non-para metric 1/V ma x dSchmidt| [l968). C + dZucca et all 1 19971) 
and STY ( San dage et al.lll979l) methods and all are in agreement 
over the considered mass range. 

Fig. [19] shows the total galaxy mass functions in redshift 
slices from z - 0.2 to z — 2. Blue circles show the results 
from the CFHTLS/WIRDS Dl field, green diamonds the D2 
(COSMOS) field, yellow triangles the D3 field and red squares 
the D4 field. Each of the four CF HTLS/WIRDS fiel ds are in 
good agreement with the results of lllbert et al.l (1201 Ol) who de- 
rived mass functions for the full COSMOS field (using the ful l 
COSMOS filter set) and with the results of iFontana eTatl {2006) 
from the GOODS-MUSIC sample. The WIRDS stellar mass 
functions shown here are presently one of the most robust mea- 
surements of the mass function at z > 1 . Taking the mass lim- 
its estimated above, the results cover a broad range of stellar 
masses, whilst the incorporation of four distinct fields allows us 
to make a robust estimate of field-to-field scatter. 

We note that, as ob s erved by a number of previous author s 
(e.g. iBaldrv et al.ll2008l: lllbert et al.ll2010t iPozzetti etalJl2010h . 
when the galaxy stellar mass function is measured over a broad 
mass range such as we have, it is not well fit by a Schechter 
function. In order to fit the observed mass functions, we there- 



fore implement the five parameter double Schechter function fit 
implemented bv lBaldrv et al.l d2008l) . This takes the form: 

where <f>(M)dM is the number density of galaxies with mass be- 
tween M and M + dM and M*, a\, a-i, <p\ and <p* 2 are the free 
parameters in the fit. We note that for the slopes, we maintain 
ci > CK2 and that at > -2, ensuring that the second term domi- 
nates at low-masses and that the luminosity density does not be- 
come divergent. The resultant double Schechter function fits are 
plotted (long-dashed curves) in Fig. [19] for each redshift/field 
combination. The double Schechter function fits consistently 
provide reliable fits to the mass function. 

In Fig. [20] we show the double Schechter function fit param- 
eters as a function of redshift for our sample. The blue, green, 
yellow and red lines in each panel show the parameters for the 
fields Dl, D2, D3 and D4 respectively, whilst the hatched re- 
gions show the uncertainties on the parameters in each case. The 
parameter values are also given in Table [7] 

Using the Schechter function fits, we calculate the galaxy 
number densities for given minimum mass cuts as a function 
of redshift. The results are shown in Fig. [2T] where we show 
the number densities for mass cuts of M > 10 9 25 M Q (blue), 
M > 10 9 75 M o (cyan), M > 10 10 25 M o (green), M > 10 la75 M G 
(yellow)andM > 1O U M (red). In each case, the solid lines give 
the mean number density across the four fields, with the shaded 
regions giving the scatter (ltx) between the four fields. Short- 
dashed lines denote results where M„,,„ is less than the mass 
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Fig. 20. Evolution of the five parameters in the double Schechter 
function fits (M* , a\, a%, 4>* and (p*) as a function of redshift. 
The hatched regions show the Icr uncertainty on the fit parame- 
ters. 



completeness limit for a given redshift. We compare our results 
to a number of other surveys denoted by the points with error 
bars. These are calculated in the same ma nner as the W IRDS 
data usin g the published Schechte r fits for lPozzetti et al.l d2007l 
squares), Marchesi ni et alj j2009i upside-down triangles) and 
lllbertetal J (1201(1 triangles) a nd from the published point s in the 
case of Fontanaetal. (2006, ci rcles) and iPozzetti et alj d2010L 
stars). The results of lllbert et ail d2010h and lPozzetti et alj (12010. 
stars) and those from the WIRDS D2 field (long dashed lines) 
are all obtained from the same field (COSMOS), although with 
slightly differing datasets. It is promising to see that these are 
all consistent with each other as shown in Fig. [21] Additionally, 
an important result from this analysis is that we see the field- 
to-field variation between the WIRDS results in the four fields, 
which is highlighted by the high values for the number density 
in the COSMOS field at redshifts of z ~ 0.3 and z ~ 0.6 - 1. 
These high densities are seen in all three D2/COSMOS measure- 
ments and relate to known large scale over-densi t ies in the field 
dFinoguenov et all 120071 iMcCracken et al J 120071: iMeneux et alj 
2009). By combining the four WIRDS fields, we see a more 
accurate picture of the stellar mass content traced over cosmic 
time. 

The WIRDS data are consistent with mass assembly down- 
sizing in which the more massive galaxy populations (i.e. M > 
1O 1O75 M ) are mostly in place by z ~ 1, whilst the numbers 
of lower mass objects are still increasing with lower redshifts. 
For instance, increases of 0.19 dex, 0.19 dex and 0.08 dex are 
seen for the mass limits of 1O 9 25 M , 10 9 ' 75 M o and 1O IO25 M 
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Fig. 21. Evolution of the galaxy number density as a function of 
redshift for the mass limits log(M/M ) > 9.0, log(M/M ) > 9.5, 
log(M/M ) > 10.0, log(M/M ) > 10.5 and log(M/M ) > 11.0. 
For each mass limit, the solid line gives the mean number density 
across the four fields, whilst the shaded region gives the scatter 
between the fields. Dashed lines give the measurements that are 
incomplete based on our mass-completeness estimates. Below 
Z ~ 1, the numbers of massive galaxies remain relatively con- 
stant, whilst the numbers of low-mass galaxies increase with de- 
creasing redshift. Above z ~ 1, the build-up of massive galax- 
ies is seen with a significant increase in galaxies of masses 
log(M/M ) > 10.5 and log(M/M ) > 1 1 .0 evident from z ~ 2 to 
z~ 1. 



respectively from z = 0.8 to z = 0.3. No increase is evident in 
the number densities between z = 0.8 to z = 0.3 in the higher 
mass samples. Comparing to data from the literature, we find 
a good consistency between the datasets, in particular for the 
lower mass cuts. At higher masses, the range of data consistently 
shows the number of massive galaxies reaching a maximum at 
z ~ 0.7 - 1.2. Although the size of our fields means we are un- 
able to probe much higher in mass than ~ 1O U M in this way, we 
note that other authors have reported this trend continuing with 
higher mass galaxy populations havi ng reached their prese nt day 
numbers at ever higher redshift (e.g. lMortlock et al]|201 ll) . 

This lack of evolution in the massive galaxy population is 
consistent with studies of brightest cluster galaxies (BCGs) in 
which little change is o bserved in BCG stellar masses to red- 
shifts of z ~ 1-1.5 (IStottetalJbOlOh . These observations 
corroborate the picture of galaxy evolution whereby a mass- 
depende nt shut-down of star f or mation is experience d by galax- 
ies (e.g. iHopkins et al1l2008bllat iHickox et alj|2009l) . The cause 
of this quenching of star-formation is uncertain, but likely can- 
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Fig. 22. The observed mass functions are reproduced as in Fig. [19] (blue circles, green diamonds, yellow triangles and red squares 
for the Dl, D2, D3 and D4 fields respectively), but now with the GALFORM results as comparison. The solid line in each redshift 
range shows the equivalent result in each redshift range based on the GALFORM simulation results, whilst the dashed line again 
shows the 0.2 < z < 0.4 WIRDS mass f unction repeated fro m the first panel. The dash-dot line in the 0.8 < z < 1.0 shows the 
original z — 1 mass function directly from lBower et alj (120061) . 



didates are that the gas fuelling star-formation is exhausted by 
heavy star-formation; that the gas is heated to temperatures be- 
yond which it may collapse to form stars due to the high virial 
temperature of the host halo; or that the gas is expelled by via 
feedback from the central AGN. 



4.3. Comparison to simulations 

We compare our measurements with predictions from the 
iBower et aD ([2006 ) implementation of the "GALFORM" semi- 
analytic simulations which are shown as the solid curves in 
Fig. [22] The results of GALF ORM have b een compared to 
previous observational data by Bow er et al. who show goo d 
fits to the observed m ass functions of iGlazebrook etal] (2004), 
iFontana et al.l d2004h and iDrorv et ail d2005l) . However we note 
that the WIRDS results provide a much stronger constraint on 
the mass function over the 1 < z < 2 redshift ran ge than these 
previous results. Some of the key features of the IBower et alJ 
GALFORM implementation are (i) a time scale for quiescent 
star formation that varies with the dynamical time of the disk and 
which therefore changes significantly with redshift, (ii) bursts of 
star formation occur due to both galaxy mergers and when disks 
become dynamically unstable, and (iii) the inclusion of both su- 
pernova and active galactic nuclei (AGN) feedback. This feed- 
back is implemented in such a way that active galactic nucleii 
(AGNs) are able to heat t he cooling fl ows in massive haloes, 
quenching star formation. IBower et al.l adopt the cosmological 



parameters of the Millennium Simulation (Sprin ge! et al.ll2005h . 
consistent with cosmic microwav e background radiation and 
large scale galaxy clustering (e.g. Sanchez ~et alJl2009h : = 
0.25, Ap = 0.75, Qi = 0.045, cr 8 = 0.9 and h = 0.73. The 
IBower et al.l model parameters were fixed with reference to a 
subset of the available observations of galaxies, mostly at lo w 
redshift. For further details we refer the reader to Bow er et all 

We note that GALFORM uses a Kennicutt IMF, whilst the 
photometric masses deri ved from the WIR DS data assume a 
Chabrier IMF. Based on Illbert etalJ (1201 Oh and iGilbank et all 
d2011l) . we multiply the GALFORM masses by a factor of 1.32 
to match the Chabrier ba sed stellar ma sses of the WIRDS pho- 
tometric catalogues. The IBower et alJ GALFORM results have 
been shown to successfully reproduce the stellar mass function 
up to z - 4.5 and the number counts of red gal axies at z < 2 
(I Almeida et al.ll2008l:lGonzalez-Perez et al.ll2009l) . 

The GALFORM mass functions as a function of redshift are 
plotted alongside the WIRDS results in Fig. [22] (solid curves). 
The da sh-dot line in the 0.8 < z < 1.0 panel shows the 
original IBower et al.l z — 1 mass functioifl Here we see that 
the GALFORM model shows some success in reproducing the 
CFHTLS WIRDS mass functions up to z ~ 2, showing good 



8 We note that the shape of the mass function predicted with the 
IBower et al.l model presented here di ffers slightly from that in the orig- 
inal paper. After the publication of IBower et alJ a bug was corrected, 
forcing slight changes in a number of parameters, affecting the shape of 
the mass function. 
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agreement at intermediate masses with both the overall num- 
bers and the position of the break in the mass function. We 
note however some deviations of the model from the observa- 
tions, with the counts being over-predicted at log(M/M ) < 10, 
whilst the model begins to generally under-predict the numbers 
of galaxies at masses of log(M/M Q ) ~ 11, especially at red- 
shifts of z > 1 • The over-prediction of galaxies with masses of 
log(M/M ) < 10 may be attributed to the form and strength of 
the super nova-feedback pres criptions in the simulation, as dis- 
cussed in lBower et al 1 (120121) . 

We also note that the models show an evolution with redshift 
in the slope of the mass function above M* , suggesting the rate 
of increase in numbers of galaxies with decreasing redshift in- 
creases as a function of stellar mass above M* . It is not clear 
that this is seen in the data where the slope at M > M* appears 
relatively consistent across the redshift slices where it is clearly 
defined. We also note that the GALFORM results are not con- 
volved with the observational errors and that this may reduce the 
tension in terms of the evolution or lack thereof of the high mass 
slopes. 

5. Conclusions 

In this paper we have presented a new deep, JHK S near- 
infrared survey using WIRCam which partially covers the four 
1 deg 2 CFHTLS deep fields. The addition of high-quality near- 
infrared data reaching K ~ 24 (AB) is an essential tool for ad- 
dressing a number of crucial outstanding questions in galaxy 
formation and evolution including, but not limited to, a pre- 
cise measure ment of luminosity and st e llar mass functions at 
z > 1 (e.g . Fonta na et al. 2004, 12006t ICirasuolo et all 120071; 
lllbert et alJl2010h : investigating the evolutionary history of the 
most massive galaxies beyond z ~ 3 (e.g. iFontana et alJ 
2006) and probing the changing relations hip between baryons 
and dark matter beyond redshift one (e.g. IFoucaud et al.ll2010t 
IHartlev et alJl2010l;(Wake et alj|201 lb . Furthermore, deep, wide, 
near-infrared data facilitates the identification of high-redshift 

I I I _| ■ " | 

groups and clusters (e.g. Finoguenov et al. 2007; Bielbv et al. 
2010[ iFinoguenov et alT 1201 Ol Henry et al] 120101; iGobat et alJ 
201 it iFassbender et al.ll201 lh a nd the first populations of galax- 
ies in the early Univers e (e.g. iHu & Cowl 3 120061: [Ouchi et al. 
l20lotlCapak etaT]l20Tlh . 

The final images cover areas of 0.49 deg 2 , 0.78 deg 2 , 
0.40 deg 2 and 0.36 deg 2 in the Dl, D2/COSMOS, D3/AEGIS 
and D4/LBQS2212-17 fields respectively, giving a total effective 
area of 2.03 deg 2 . Measured full-width-half maxima in all fields 
are « 0.8" or better (except for the D2/COSMOS WFCAM J- 
band image, which has a FWHM ~ 1.0". For the Dl, D2 and D4 
CFHTLS fields, the depths in all three bands are consistently bet- 
ter than magAB ~ 24.5 (50% point-source completeness), whilst 
the D2/COSMOS field is slightly shallower reaching J = 23.4, 
H = 24. 1 and K s = 24.0. We have presented combined counts 
for our 2.03 deg 2 survey in all three bands. These are consistent 
with previous works, whilst our measurements are generally ei- 
ther made over a larger area or deeper than previous measure- 
ments. 

WIRDS data have been combined with CFHTLS deep op- 
tical data in each of our four fields to create high quality 8- 
band photometric catalogues. From the full 8-band ugrizJHK s 
data-set we have constructed two sets of combined catalogues, 
the first catalogue set made using x 2 gfi detection images and 
the second set constructed using K s images as the detection 
source. All WIRDS images and catalogues are publicly avail- 
able to download from the CADC. Photometric redshifts have 



been determined based on both of these catalogue sets using the 
Le PHARE software. In addition we have produced a catalogue 
of galaxy properties based on the SED fitting using Le PHARE 
using the^ 2 gri based catalogue, giving estimates of galaxy stel- 
lar masses. 

We have presented an analysis of the success of the BzK 
selection based on the CFHTLS/WIRDS data compared to full 
photometric SED analyses for redshift and galaxy type determi- 
nation. Our analysis has shown that the pBzK selection based on 
the CFHTLS/WIRDS filters successfully selects * 52% of the 
passive galaxy population at 1 .4 < z < 2 when compared to the 
photometric SED-fitted catalogue. In addition, the selection se- 
lects a significant number of 1 < z < 1.4 galaxies. The sBzK 
selection successfully identifies 76% of the 1.4 < z < 2 star- 
forming galaxy population identified by our SED fitting, with a 
much stronger cut at z — 1 .4 than seen for the pBzK selection. 

We have presented estimated mass limits for our /-selected 
photometric catalogue and presented the total galaxy mass func- 
tion as a function of redshift over the range 0.2 < z < 2. 
These results provide one of the most robust measurements of 
the galaxy stellar mass function at z > 1 presently available. In 
particular, we have shown that the mass functions are consistent 
between each of the four individual fields covered in this sur- 
vey, providing a strong constraint on the effect of cosmic vari- 
ance on our measurements. The results are consistent with the 
current best determinations of the stellar mass function over our 
redshift range. We have shown that the mass functions are well 
fit by double Schechter function fits, whilst noting that single 
Schechter functions do not provide good fits across the broad 
mass ranges covered. The evolution of the Schechter function 
parameters have been presented as a function of redshift for each 
of the four fields, noting in particular that M* shows a grad- 
ual decline from z = to z = 2. From the fits to the mass 
functions, we have calculated the number densities of galaxies 
above a set of mass limits as a function of redshift. The results 
are consistent between the for fields and show the build-up of 
mass as a function of redshift since z = 2. We find that since 
z ~ 0.8 - 1, there has been little change in the numbers of mas- 
sive galaxies (M > 10 10 75 ), whilst the formation of lower mass 
galaxies has been ongoing in the same epoc h. This is consistent 
with th e previous fi nding s of for example IPozzetti et alJ (2007, 
1201(1 ; lllbert et all (12010) . From z ~ 2 to z ~ 1, the WIRDS 
data clearly show the significant increase in numbers of galax- 
ies of all masses probed by the survey as the Universe under- 
goes the phase of peak star-forming activity dLillv et alj fl996: 
lMadauetalJll996l) . 

Finally, we have compared our results with the predictions 
of the semi-analytical galaxy formation model GALFORM and 
find that the simulations provide a relatively successful fit to 
the observed mass functions at intermediate masses of 10 < 
log(M/M ) < 1 1 . However, as is common with semi-analytical 
predictions of the mass function, the GALFORM results under- 
predict the mass function at low masses (i.e. log(M/M Q ) < 10), 
whilst the fit as a whole degrades beyond redshifts of z ~ 1.2. 
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Table 7. Double Schechter function fit parameters (and associated lcr uncertainties) for each field as a function of redshift. 
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